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Abstract 

We carry out a systematic construction of the coarse-grained dynamical equation of motion for the 
orientational order parameter for a two-dimensional active nematic, that is a nonequilibrium steady 
state with uniaxial, apolar orientational order. Using the dynamical renormalization group, we show 
that the leading nonlinearities in this equation are marginally irrelevant. We discover a special limit 
of parameters in which the equation of motion for the angle field of bears a close relation to the 2d 
stochastic Burgers equation. We find nevertheless that, unlike for the Burgers problem, the nonlinear- 
ity is marginally irrelevant even in this special limit, as a result of of a hidden fluctuation-dissipation 
relation. 2d active nematics therefore have quasi-long-range order, just like their equilibrium counter- 
parts 



*Electronic address: smishr02@syr.edu 
Electronic address: |aditi@physics.iitni.ac.in| 
^ Electronic address: |sriram@physics.iisc.ernet"Iii| 



1 



Active nematics are the simplest example of spontaneously broken rotation- 

invariance in a nonequilibrium system. Analytical studies of their statistical properties 
have mainly been confined to a linearized approximation whose predictions of anoma- 
lous density fluctuations have largely been confirmed in experiments 3| and numerical 
simulations j^. Within the theory of pj] the density fluctuations were driven by the 
broken-symmetry modes associated with orientational order. In this paper we ignore den- 
sity fluctuations and focus on the effect of the broken-symmetry modes on the strength 
of orientational order. We ask: can a noisy two-dimensional system of active particles 
display long-range nematic order? 

Let us see why this question is worth asking. It is well known that at thermal equilibrium, 
in two space dimensions, neither XY models nor nematic liquid crystals can have long- 
range order. Instead of a true ordered phase, these systems have a critical low-temperature 
state in which the fluctuation-averaged order parameter vanishes in the thermodynamic 



limit at al 

distance 



nonzero temperatures, but order-parameter correlations decay as a power of 



7(1 . The simplestgeneralization of the 2d XY model to a nonequilibrium 
steady state is the Vicsek model 8|] of flocks in two dimensions, in which the local velocity 
of the flock is the XY order-parameter field. Toner and Tu |9| showed that the resulting 



advection of the order-parameter field by its own fluctuations jo] stabilizes long-range 
order even in two dimensions. Technically, the mechanism amounted to a singular renor- 
malization of the XY stiffness by nonlinearities of a type not permitted in the equilibrium 
XY model. The ordered state of a Vicsek flock can be thought of as a collection of arrows 
all pointing on average in the same direction; this is known as polar order. One can imag- 
ine a different ordered state, in which the axes of the arrows are on average parallel to 
an arbitrarily chosen spatial direction, call it n, but the arrows point indifferently along 
-|-n and — n or, equivalently, one could simply lop the heads off the arrows. The resulting 
state is apolar, and has purely nematic order. The Vicsek flock moves on average in the 
n direction, while a nonequilibrium steady state with nematic order - an active nematic 
- cannot tell forward from back, and so does not drift on the average. The nature of 
order in such active nematics is the subject of our study. Our main concern is whether 
the interplay of nonlinearity and fluctuations stiffens the order-parameter fluctuations in 
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active nematics as it does 9| in polar ordered phases, leading to true long-range order in 
two dimensions. 

Here are our main results, (i) We elucidate the route to the equation of motion for the 
nematic orientational order parameter, taking care to distinguish the constraints intro- 
duced purely by rotation invariance, and hence applicable to both active and equilibrium 
systems, from those arising specifically in the thermal-equilibrium limit, (ii) We show 
that the two quadratic nonlinearities in the equation of motion have independent coeffi- 
cients, unlike in the equilibrium case where they are determined by a single parameter. 
In both equilibrium and active nematics power-counting shows that the nonlinearities are 
marginal, but such analysis cannot distinguish marginally relevant from marginally irrel- 
evant, (iii) In a certain limit of parameter values, our equation of motion can be mapped 
to the noisy two-dimensional Burgers [l^ and KPZ 11| equations, but with a veloc- 
ity field V satisfying the peculiar condition d^v^ — d^Vz = 0, which is neither solenoidal 
nor irrotational. (iv) The similarity to the Burgers problem ends there: our dynamical 
renormalization-group treatment shows that the nonlinearities are marginally irrelevant 
in our theory, in the Burgers limit as well as in general. Active nematics thus have only 
quasi-long-range order. Although disappointing if one is looking for novelty in nonequi- 
librium systems, this negative result reinforces the findings of a numerical study of an 
apolar generalization of the Vicsek model. 

This paper is organized as follows. In sectionUwe construct the coarse-grained equations of 
motion for the nematic order parameter, highlighting the differences between equilibrium 
and active systems. In section II Bl we examine the relation of our equations of motion 
to the Burgers and KPZ equations, in a special high-symmetry limit. In section [TTl we 
outline the dynamic renormalization group (DRG) treatment with which we extract the 
long-time, long-wavelength properties of correlation functions in our system. Further 
calculational details are relegated to the Appendix. The paper closes in section IIIII with 
a discussion of possible future directions. 
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I. EQUATION OF MOTION 



We now construct the equations of motion for an active nematic. Since we are considering 
a system that can undergo apolar orientational ordering, one of the slow variables for a 
coarse-grained description of the dynamics is the traceless symmetric second-rank tensor 



nematic order parameter Q [12|]. The magnitude of Q is slow upon approach to the 
ordering transition, and the fluctuations of its principal axis are the broken-symmetry 
modes of the ordered phase. If the system were isolated, mass and momentum would be 
conserved within the system and the corresponding densities p and J = pv, v being the 



velocity field, would be slow variables as well [13|. However, we will consider a system 
adsorbed on a solid surface which acts as a momentum sink, thus turning J or v into a 
fast variable, and allow deposition and evaporation i.e., birth and death 15], thus 



rendering p fast as well. We will start from a complete dynamical description, eliminate 
the fast p and J, and obtain the dynamics of Q alone. 

For a system where particles can enter and leave the system in the bulk, the density obeys 

|^ = -7P + /?-V-J + /,. (1) 

The third term on the right of ([T]) contains the number-conserving motion of particles on 
the substrate. The random adsorption and desorption of discrete particles has two effects. 
In the mean, conditioned on a given local density p(r,t), it leads to the 7 and (3 terms. 
Fluctuations about this average effect lead to the nonconserving spatiotemporally white 
noise fp. A steady, spatially uniform state has mean density po = /5/7- Newton's second 
law for the momentum density mJ reads 

m— = -Fv + - V ■ a (2) 

The first term on the right hand side of ([2D is friction due to the substrate, with a kinetic 
coefficient F. The random agitation of the particles as a result of thermal motions, 
biochemical stochasticity, or dynamical chaos is modelled in the simplest possible manner 
by the spatiotemporally white Gaussian noise f/j. This noise is nonconserving, i.e., its 
strength is nonvanishing at zero wavenumber, since the dynamics is not momentum- 
conserving. The last term contains all effects arising from interactions of the particles 



with each other, and thus takes the momentum-conserving form of the divergence of a 
stress tensor a. In principle a contains stresses coming from the free-energy functional for 
Q (see below) [l^ These, however, are readily seen [l| to be irrelevant at large lengthscales 
compared to the contribution cr" = WipQ coming from the active nature of the particles 



17|. 



The equation of motion for the orientational order parameter Q including coupling to the 



velocity field 
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19| is 



dQ 

_ + V ■ VQ = rC + {aoK + aiK ■ Q)st + n Q-Q n (3) 



where k = [Vv+(Vv)"^]/2 and Q = [Vv — (Vv)"^]/2 are the shear rate and vorticity tensor 
respectively, F is a kinetic coefficient [2^ , and the parameters and ai characterise the 
coupling of orientation to flow. The molecular field G = —6F/6Q is obtained from an 
extended Landau-de Gennes free energy 

F = J d'xi^TrQ' + ^{TrOr + |(V.g.O' 

+ KQijViQkiVjQki + CQi.ViVjp] + $[p] (4) 



where we have left out terms cubic in Q as these vanish [1^ in dimension d = 2. The 
density p enters F through the functional $, the quadrupolar coupling term with coeffi- 
cient C, and the p-dependence of parameters in /. On timescales much larger than I/7 
and m/r, the density and momentum equations ([T]) and ([2]) become constitutive rela- 
tions determining p and J in terms of the slow field Q. Eq. ([T] tells us we can replace p 
everywhere by po to leading order in gradients, and becomes 

--■ "^^°-V-Q (5) 



r 

apart from noise terms. The molecular field G in (I3l) contains a term of the form QVVQ, 
and one of the form VQVQ, whose coefficients will be related as both terms arise as 
variational derivatives of the single K term in F [Eq. (jlj). Replacing v by its expression 
([5]) in Eq. (jSj) will give rise to additional terms of that form, controlled by the activity 
parameter wi. As a result, the QVVQ and VQVQ terms in the effective equation of 
motion for Q cannot be combined into the variational derivative of a scalar functional. 



and will have two independent coefficients. We will explore below the consequences of 
the existence of two independent nonlinear couplings. In space dimension d = 2 the 
order-parameter tensor has the simple form 

S I cos 20 sin 20 \ 
^ \ sin20 -cos 20 / 

where the scalar order parameter S measures the magnitude of nematic order and is the 
angle from a reference direction. Let us work in the nematic phase, where we can take 
S = constant and define = along axis of mean macroscopic orientation. Eq. for 
small becomes 

^ = _Y"\dAdJ), (7) 

neither a gradient nor a curl, f being a constant determined by those in ([1]) - ([6]). Sub- 
stituting V in ([3]) by its expression ([7]), writing Q in terms of as in ([H]), treating S as 
constant, and including noise terms, we obtain 
(90 

— = A.dle + A^dle + XidM.e + \2ed^d,e + fe (8) 



to order 0^, where the additive 



211 ] non-conserving Gaussian white noise satisfies 



< /e(r, t)fe{v\i) >= 2D,6{v - v)5{t - t ) (9) 

with a noise strength Dq. All the coefficients in ^ and ([9]) are related to those in ([T]) - (I3j), 
the corresponding noise strengths, and the scalar order parameter 5*. As a consequence 
of rotation invariance, i.e., the fact that the underlying equation of motion in terms of Q 
has a frame-independent form, we find 

2{Ai-A2) = \2. (10) 

It is therefore convenient to re-express them as 

= + A2/4; = Ao - A2/4 (11) 

Without the detailed derivation above, it would have been hard to guess the form of the 
equations of motion and the constraints on the parameters. Note that Ai and A2 are in 
general independent, as we argued above. We will comment below on the relation they 
satisfy in the special case of an equilibrium nematic. Eqs. (El) and ( fTOl) can also be 



obtained from a microscopic model of coUisional dynamics of apolar particles 22 1 . 



A. Equilibrium limit 



The energy cost of elastic deformations and, hence, the thermal equilibrium statistics 
of configurations, of a two-dimensional nematic are governed by the Frank free energy 
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23 



24| 



H = l [^(V-nf + ^(Vxn)V^ (12) 

a functional of the director field n = (cos^, sin 6'), with splay and bend elastic moduli Ki 
and K^. To cubic order in ^(r) 

H/ksT = ^ [ d2r[[dj{r)]^ + (1 + A)[dAr)]' 



2 

-2Ae{T)[dJ{r)d,9{r)]] (13) 

where A^, = K^/ksT and A = The purely relaxational dynamics of the angle 

field 6, at thermal equilibrium consistent with ( fT3l) . reads 

— = A^dle + (1 + A)A-idle + Xid^ed^e 
+ x^ed^d^e + /, (14) 

where {fe{v^t)fg{0^0)) = 26{r)6(t), and a kinetic coefficient has been absorbed into a 
time-rescaling. The nonlinearities in iHM have the same form as in ([8]), but the couplings 
are not independent: 2Ai = A2 = —2^3, since both come from the same anharmonic 
term in the free energy f|T3|) . In addition, the nonlinearity is connected to the diffusion 
anisotropy: 2[A3 — (1 + A)A3] = A2 as required by rotation invariance. Eq. (fT4l) is simply 
the limit 2Ai = A2 of (E]). 



A static renormalization-group treatment of the 2d equilibrium nematic [2J] with Hamil- 
tonian ( |T3i) showed that A was marginally OTelevant, and that the large-scale behaviour 
of the system was governed by a fixed point with A = 0, i.e., a single, finite Frank con- 
stant for both splay and bend. The dynamics of the active nematic does not correspond 
to downhill motion with respect to a free-energy functional, and the two nonlinear terms 
thus have independent coefficients. Their (marginal) relevance or otherwise must be es- 
tablished by a dynamic renormalization-group study of the equation of motion (IHl), which 
we present in section HIl 
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B. Burgers equation 



The structure of ([8]) in a certain special limit merits some attention. If we switch off 
the A2-nonlinearity, equation has a higher symmetry than in general, viz., under 
6^6 + constt without a corresponding transformation of the coordinates. In addition, 
it is invariant under x ^ z, which allows us to rescale the equations so that the diffusion 
of 9 is isotropic: 

— = AV'e + XdJdJ + fe (15) 

with a spatiotemporally white noise fg as in ([9]). This equation for X ^ cannot corre- 
spond to an equilibrium system, because the sole surviving nonlinear term XdxOdzO cannot 



be written as 6A/66(x.) for any scalar functional A[6] 25| Note the similarity of ( |T5l) to the 



KPZ equation ll| for the height field of a driven interface. Extending the analogy, it is 



easy to see that the velocity field v = {d^O, d^O) as in ([7j) obeys the Burgers-like equation 
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^ = AV\ + X{vV)y + U (16) 
ot 

with a conserving noise fv = {dzfe, dxfe)- The curl-free condition of a traditional Burgers 
velocity field is replaced in our case by dxVx — dzVz = 0, which amounts to equal extension 
rates along x and z. In the 2d randomly-forced Burgers-KPZ problem, the nonlinearity 



is known 
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ll[ | to be marginally relevant, so that the large-scale long-time behaviour is 
governed by a strong-coupling fixed point inaccessible to a perturbative RG. It is natural 
to ask what happens in the seemingly similar problem at hand. 



1. Galilean invariance 

Eqns. (fTSl) and (fT6l) are invariant under the infinitesimal Galilean boost 

x ^ x — ut (17) 

e^e + u-x (18) 
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or equivalently 



u 



(19) 



where 

u = (m^,m^.) (20) 

inverts the vector components of u. By analogy to the results of [l^ and H] this invari- 
ance implies that the nonlinear-coupling A does not renormalise in this special limit. 



II. RENORMALIZATION GROUP THEORY 



In this section we outline our one-loop dynamic renormalization group (DRG) analysis 
of the large-scale, long-time behaviour of Eq. ([8]). Our treatment is general, allowing 
for two independent coupling strengths Ai, A2, but we will examine the A2 ^ limit of 
section II Bl as well. We present only the key steps of the calculation, relegating details to 
the Appendices. 



The momentum-shell dynamical renormalization group (DRG) 10, l26|, Um, l28|] consists 
of two steps. Consider a system with physical fields described by Fourier modes with 
wavevector q with < g = |q| < A, the ultraviolet (UV) cutoff. First: eliminate modes 
with Ae~' < g < A, by solving for them in terms of those in < g < Ae"' and the noise, 
and average over that part of the noise whose wavenumber lies in [Ae~', A). Second: 
rescale space, time, and dynamical variables to restore the cutoff A and to preserve the 
form of the equations of motion to the extent possible. The result is an equation of 
motion in which the parameters have changed from their initial values, call them {A'o}, 
to /-dependent values {K{1)}. Now, correlation functions at small wavenumber can be 
calculated either from the original equations of motion or from those obtained after the 
above two steps. This key observation leads to a homogeneity relation between correlation 
functions 

C{ci,iu;{Ko}) = e-f'Ciqe\Lue^';{Kil)}). (21) 
that can be used to calculate long-wavelength correlations with particular ease if the 
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couplings flow to a small fixed-point value {K{oo)} under iteration of the above transfor- 
mation. Let us carry out this process for our model, Eq. ([8]). 



We insert the decomposition [29| 9{r, t) = f^^^ ^ 9{q, u) exp (iq ■ r — iut) into ([HD to ob- 
tain the 6 equation in Fourier space: 

e{q, u) = Goiq, u;)/e(q, uj) - Go{q, lu) 

I M(k,q-k)^(k,^])^(q-k,u;-fi) (22) 

JkO. 

where 

Go(q,^) = [-IUJ + A^ql + A2g,2]-i (23) 

is the bare propagator, 

M(k, q - k) = ^[K{q, - h) + hiqx - h)] 

+ y [Kh + {Qx - kx) (g. - k,)] (24) 

the bare vertex, and the Fourier transform fg[q,Lj) of the Gaussian spatiotemporally 
white noise in (IHl) has autocorrelation 

{f9{q,u;)fe{q,uj')) = 2Do{2nf^'5{q + q)6{uj + uj') (25) 



> > X 

Go(q,co) G(q.co) 4(q,co) ^ 

(q, 03 ) (q, CO ) (q, CO ) ^ 

> O < 'o. 

(a)Go(q'Co) Do G„(-q,-co) -^y 

FIG. 1: (a) Definition of symbols, (b) Diagram for full non-linear equation (j22|) in Fom'ier space. 
The left hand side of the pictorial equation is the full solution to 6'(q, u) = G'(q, a;)/e(q, w), where 
G(q, cj) is the full propagator. The first part on the right hand side is the zeroth order solution 
to (p2|) 6{q,uj) = Go(q, u;)/6i(q, w) and the second term is the contribution of the nonlinearity. 



Eq. (122!) can be represented graphically as in Fig. [H A perturbative approach to solving 
( !22|) generates corrections that can be expressed in terms of Feynman graphs of three 
types - propagator, noise strength and nonlinearities - given in Fig. O 
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(a) 




(k, 



. . ^ I \. i \. — @ — - — Q — 2 

(q,(0) (q,(0)(q-k,0]-(l) (b) 



-qrco) 

(q-k,co-fl) (k-q,fl-co) 



(c) 



+ 4- 



+ 4 { w+4 



FIG. 2: (a) Graph for propagator G(q, w). The left hand side with a double line is the full 
propagator, the first term on the right hand side is the zeroth order and the second term is 
the one-loop correction, (b) Graph for force density Z?(q,tj) defined by (|25p. The second 
term on the right hand side is the one-loop correction, (c) Graph for the three-point vertex 
function. The structure with three legs with one incoming and two outgoing is the vertex 
— ^2tt)'^+i I -^(^j q ~ k). The three graphs are Fq, F^ and Fc. 



A. Propagator calculation 



The effective propagator G(q, cu) [defined by 9{q,uj) = G{q,uj)fe{q,uj)] is given pertur- 
batively in Fig. [2](a). The averaging over the noise is performed using fl25|) . The one- loop 
correction to the propagator is 

^(q, u) = Go(q, ^) + 4G'2((q, uj) x 2Do 

I M(k,q-k)Af(-k,q)Go(k,fi) 

Go(-k,-f^)Go(q-k,cu-fi) (26) 

or 

^-^(q, u) = Go ^(q, OJ) - S(q, u) (27) 

with a self-energy 

S(q, u}) = Ax2Dq I M(k, q - k) M (-k, q) 

Jkn 

G'o(k,f})Go(-k,-fi)Go(q-k,cu-f]) (28) 
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where the combinatorial factor of four represents possible noise contractions leading to 
Fig [2] (a). A few steps of calculation of the integrals are performed in Appendix |X1 For 
small wavenumber q and for — > 0, the result of integrating out a shell between Ae~' 
and A in q space is the self-energy. 



S(q,0) 



+ 



^ _ (j2(Ai, A2) 

8 

G3{Xi,X2)AiA2 



(29) 



where 



G^2(Ai, A2) — (2A^ + A2 — 3A1A2 
CslAi, A2) = (A2 — A1A2) 



The dimensionless quantities Ai and A2 are defined by 

AjAj = — — / /9 ' ^'-^ ~ -'-'^ 



(30) 



(31) 



(^1^2)3/2' 

When we implement the dynamical renormalization group, terms of order and of order 
1 are generated though the self-energy. Terms of order will give corrections to the 
diffusion constants (Ai, A2). What about the terms [3^ of order 1, which also arise in the 
analysis of Pelcovits et al. 24]? As in 24], we proceed by first ignoring the terms of order 
1, whose coefficient is proportional to one nonlinear coupling A2, and then, post facto, 
realise they too are (marginally) irrelevant because A2 itself is found to be marginally 
irrelevant. Proceeding in this manner we find 

G-i(q,0)=Go^(q,0)-S(q,0) 

~ A^ql + A2ql 

= A^ql + A^ql 
, G2(Ai,A2)(Ag2 + A2g2)/ 



That is, 



ii = Ai[l 
^2 = ^2(1 



4 X Svr 



G'2(Ai, A2)/-, 
4 X Svr ^ 
G'2(Ai, A2)/i 



(32) 



(33) 



4 X Svr ' 

These are the intermediate (one-loop graphical) corrections for anisotropic diffusion con- 
stants. 
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B. Vertex calculation 



From the full equation (1221) and (Fig[T]), the diagrams contributing to the vertex correction 
are shown in (Fig[2](b)). There will be three types of diagrams, all with multiplicity 4, 
denoted by F^, F^ and F^. The details of the calculation are given in Appendix [Bl The 
full vertex is defined as a combination of Ai and A2 equation (!24l) . We study how this 
vertex evolves under the DRG and at the end of the calculation we can separate terms 
corresponding to Ai and A2. From (Fig[2](b)), expression for 

F„(q,ki) = 4x 2Do / M(k,q-k) 

Jkn 

xM(| + ki,k-|-ki) 
xM(|-ki,-k+| + ki) 

2 

X Go{k,n)Go{q-K^ (34) 

The integral as usual is over Ae^' < g < A. Similarly one can get expressions for Ffe(q, ki) 
and Fc(q, ki). Hence, adding contributions to all diagrams for the vertex, Fa(q, ki) + 
F6(q, ki) + Fc(q, ki) we can get the graphical corrections to the couplings Ai and A2. 
After a calculation as in Appendix [Bl the graphical corrections to Ai and A2 are 

Fi(Ai,A2)/, 



Ai = Ai[l 



2 X Stt 



i^i(Ai,A2),F2(Ai,A2) defined by, 

i^i(Ai, A2) = -2A1A2 + 3X1 + ^2Ai 

F2(Ai,A2) = -4A2A1 + 6A^ (36) 

Note from ( l36i) that Fi(Ai,A2) = 0, if A2 is zero. This says that there is no graphical 
correction to Ai if A2 is zero. This is a result of the Galilean invariance in this limit, as 
pointed out in section II B 1[ 
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C. Noise strength renormalization 



An effective noise strength D can be defined by 



r (q, uj)d{ci, uj)) = 2DG(q, ^)G'(-q, -uo) 



(37) 



This quantity is calculated perturbatively by the series shown in (Fig[2]^c)). To one- loop 
order 

2D = 2Do + 2{2DoY [ M(k, q - k)M(-k, k - q) 

Jkn 



X 



Go(q — k, UJ — Q) 



(38) 



The integral in equation 0381) is performed in Appendix O After doing the integrals, the 
graphical correction to Dq is 

, (A2-Ai)2/- 



D = Dn 



2 X 87r 



(39) 



1. The detailed balance limit 



From equations ([33]) and ([39]), for Aa = {Ai = A2 = A, G2(Ai,A2) = (2A? + A2-3AiA2) = 
2A^ and (A2 — Ai)^ = A^), i.e., A and D have the same graphical corrections. This suggests 
that detailed balance should obtain in the limit A2 = 0. To discover this detailed balance 
let us write the Fokker-Planck equation 3l| for the probability distribution functional 
P[e,t] of the e-field: 

d 



dP ^ d 



+ 



^5];M(1, m)^l^^5q,l 
* l,m 

We guess that a Gaussian probability distribution function 



P = 0. 



(40) 



(41) 
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is a steady solution to equation (HOj) . M(l, m) = {Ixfriy + rrixly), is a normalization 
factor and the two-point function < 6^0 >= {Dq/ A)q~'^ . If this is so, the last term on 
the right of equation (l4Up should vanish if Pgt from equation (HTj) is inserted for P. Let 
us check this: 



q,l,m ^ 



dPo 



^M(l,m)^,^^(5q,i+, 

q,l,m ^ 

-PO^ J2 m)^l^m^-q5q,l+m (42) 

q,l,m 



Using the symmetry — q ^ 1 ^ m in fH2]) we get 

m 

- ^[M{\, m)(l + m)2 + l^M{ — m, 1 + m 



1 

3 

m^M(-l,l + m)]eie^^_i_„ (43) 



l,m 
.2 ; 



The summation inside the square bracket in (143!) is zero. This means that for A2 = the 
Gaussian defined in (HTl) . is a steady solution of the FP equation PUI) . consistent with 
the detailed balance noted after equation ( 1391) in this limit. In particular, we can already 
conclude that there is no singular renormalization of the stiffnesses in the Burgers-like 
limit of the model, as the equal-time correlators of 6 can be obtained directly from the 
Gaussian probability distribution function (14T|) . 
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D. Full RG Analysis 



We now return to the general case Ai, A2 nonzero. Substituting results from (1331) . (135!) 
and fl5^ to ([22D, gives the intermediate equation for 6''^(q, cj) (without rescaling) 

-Giiq,L0) [ Mi{k,q-k) 



JkQ 

X e<{k,n)e<{q-'k,uj -n), (44) 

where the propagator at this intermediate stage is 

Gi{q,u;) = i-tu + A^ql + A^ql)-', (45) 

with Ai and A2 given by fl33l) and < |q| < Ae~', unlike the original equation, which is 
defined on the large range < |q| < A. 

Next rescale variables to preserve the form of the original equation: 

q = ge'; lu' = cje"^'^; 

e<{q,u) = m0'i^.^')- (46) 

Thus the new variable q' is defined on the same interval < |q'| < A as the wave- vector 
q in the original equation. In terms of the new variables, the intermediate equation for 
9'{q_',uj') is 

^'(q',^') = G'(/)(q',^')/e(q'7^') 

-G(/)(q>') / M(/)(k',q'-k') 
Jk'n' 



where 



with 



X ^ (k,a)^ (q -k,cu -a), (47) 

G(/)(q>') = Ha;+Ai(0g;^+^2(09:r' (48) 

Ai(/) = iie°«-2'; A^il) = iae"^'^-''; (49) 

M,uj') = f<{q,u;)e'^^'^r\l) (50) 
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M(0(k',q'-k' 



Ai(0 



+ 



X 



(51) 



k'X + (Qx - K)iQz - K) 
where Ai(/) and A2(/) are rescaled nonlinearities given by 

Ai(/) = Aie(/)e-('^+2)^. ^ ^^^^^^^_(,+2); ^52) 

The correlation function characterising the force /g(q',u;'), given by expression (|50l) . can 
be constructed using definition ( 125|) and the new set of variables ( l46l) 

< /;(q,.;)/;(q>') >=2i;(0 (27r)2+^5(q + q') 

(5(c<j + c<j') (53) 

with 

D{1) = £)e(=^"(')+'^')r'(0 (54) 
where d = 2 and all tilde variables correspond to the graphically corrected quantities in 
(133|1 . fl35|) and fl39l) . Substituting for the expressions for all tilde variables 

Ai(/) = Ai[l + ^^^%M]e"{0-2/, 

A2(o = ^2[i + ^^^|^r(')-^', 

D(/) = Z^[l + i(^^]e-(0-^-r^(/). 



(55) 



E. Recursion relation 



Here we calculate the recursion relation for all five parameters. From (I55p . the constraint 

of rotational invariance 2{Ai — A2) = A2 requires 

I n\ I , ^^2(Ai,A2) /-F2((Ai, A2) \ . . 

m = eMail) + 21) [1 + ,^3^ + 4^8vr ) ^''^ 
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where the functions 6*2(^1, A2) and F2(Ai, A2) are aheady defined in (130|) and (136!) . With 
this choice of substituting in fl55l) . recursion relations for all five variables given by, 



dAi 
"df 
dA 
dl 
dXi 
'df 



Ai[-2 + z{l) + 



: A2[-2 
Ai[-2- 



-zil) 
z{l)- 



G'2(Ai, A2) 

4 X Stt 
G'2(Ai, A2) 



4 X Stt 
F2(Ai, A2) 
4 X Stt 



+ 



G'2(Ai, A2) 

4 X Stt 



i^i(Ai,A2) 
4 X Stt 



dX2 
'df 



\2[-2 + z{l) + 



^ = D[-2 + z(l) + 
dl ^ ^ ' 2 X Stt 

_ -p2(Ai, A2) . 

2 X Stt ^' 



G'2(Ai, A2) n 
4x87r ^' 
(A2-Ai)2 G'2(Ai,A2) 



2 X Stt 



(57) 



where z{l) is defined by a{l) = z{l')dl' , and the dimensionless variables Ai and A2 were 
defined in ( !3T|) . The functions G2, Fi and F2 are already defined in (130|) and (1361) . In 
these recursion relations the function z{l) is unknown at this point. It will drop out in 
the recursion relation for the dimensionless variables, Ai and A2, for which the recursion 
relations are 



dXi 

~dr 

d\2 
IT 



Ai 



A. 



(A2 



Ai)^ 3 G2(Ai,A2) 
2 4 X Stt 



i^r(Ai,A2) 
4 X Stt 



4 X Stt 

(A2-Ai)=^ _ 3 G2(Ai,A2) _ F2(Ai,A2) 
2 4 X Stt 



4 X Stt 2 4 x Stt 4 x Stt 

Equations (l58l) and ( l59l) are coupled nonlinear equations for Ai and A2. 
In the special, high-symmetry case A2 = 0, from ( l30l) . G2(Ai, A2) = 2Af and -Fi(Ai, A2; 
Then the dimensionless coupling A^(/) = Xl{l)D{l)/A^^'^{l) obeys 



dXi 
'df 



Ai 



-2 + 2 



m 

2 X Stt 



m 

2 X Stt 



(58) 
(59) 

= 0. 

(60) 



which tells us Ai is marginally irrelevant. By contrast, for the Burgers equation in 2-d, 
the nonlinearity is marginally relevant. This is surprising, given the similarities of the two 
models in the limit A2 = 0. A second special case is A2 = 2Ai, when the problem reduces 
to an equilibrium problem, as remarked in section II A[ At this particular choice of Ai 



and A2, G2(Ai,A2) = 0, (A2 - Ai)^ = A? = f , Fi(Ai,A2) = 16A? and F2(Ai,A2 



4A^. 



Substituting these expressions for all functions in (158|) and (l59l) . the flow equations for 
the equilibrium limit are 

rf/ 4 X Svr ^' c// 4 X 4 X Svr 2 ^ ' 

We can draw the flow-diagram in (Ai, A2) plane. (Fig[3]) shows that for three special cases, 
A2 = 0, A2 = 2Ai and A2 = Ai flow is towards zero. For other points also flow is towards 
zero. This means (0, 0) is the only fixed point and it is stable. We have checked this 
numerically as well. 

Since the nonlinearities are marginally irrelevant the effective stiffness Ai and A2 become 
equal at large scales, and are nonsingular. Therefore < |6'gp >~ for small q, i.e. the 
renormalized theory still has only quasi long-ranged order. 



III. CONCLUSION AND DISCUSSION 



In this paper we have provided a systematic analysis of the large-scale, long-time behaviour 
of the stochastic nonlinear partial differential equation for the angle field of an active 
nematic on a 2-dimensional substrate. We constructed the general equation of motion 
for the order parameter, starting from a description that included the velocity, density 
as well. We then reduced the model to focus on the director or small-angle fluctuations 
about an ordered active nematic, and studied the evolution of the parameters therein 
under the dynamic renormalization group [lO, 11, 26|. The equation has five parameters. 



Ai and A2 which are director diffusivities for two directions, the nonlinear couplings Ai 
and A2 and -Do the noise strength. Two special cases are of interest: A2 = 2Ai, for which 
the dynamics is that of an equilibrium two-dimensional nematic where static properties 



are shown to agree with [2J]. The second case is A2 = 0, for which the equation can be 
mapped to a Burgers equation, for a velocity field v given in ([7]), with d^Vx — dzVz = 0. 
Despite this resemblance the dimensionless nonlinear coupling parameter A^ = i s 



io|) 



found to be marginally irrelevant, whereas for the Burgers equation in d = 2 (see 
the nonlinearity was marginally relevant. Interestingly in this limit the diffusion constant 
and noise strength renormalize the same way, implying the system has a hidden detailed 
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balance, which we exposed via a Fokker-Planck analysis. The complete one-loop recursion 
relation for the five parameters constrained only by rotational-invariance show that the 
nonlinearities are always marginally irrelevant. 
In Appendix |D] we present the equation of motion for the angle field starting from 




FIG. 3: RG flow diagram in the phase plane of dimensionless nonlinear couplings Ai and A2 
defined in (j3ip . The solid line represents line A2 = 2Ai (equilibrium limit), the dot dashed line 
represents Ai = A2 and dashed line represents A2 = (limit when equation is similar to Burgers 
equation). For these three cases, it is particularly easy to show analytically that the flow is 
inward (i.e. nonlinearities are marginally irrelevant). In fact for all Ai, A2 the flow is towards 
(0, 0). 

a velocity field which satisfies incompressibility. This provides another, inequivalent, 
situation in which the density is fast and can therefore be suitably eliminated. The 
procedure leads to a slightly different equation from ([8]) or (1221) with nonlocality due to 
transverse projectors. We have not analysed the properties of the incompressible version. 
Our results, despite the neglect of the density, are consistent with the numerical findings 
of [2I, that active nematic order in (i = 2 is quasi long-range. A complete treatment of the 
coupled behaviour of angle and density correlators in steady state, beyond the linearized 
analysis of as well as a study of the incompressible model, are left for future work. 
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APPENDIX A: PROPAGATOR RENORMALIZATION 



We start from the symmetrised version of (12^ (by substituting k = | + k and f2 = ^ + f2) 

E{q,u) = 4x2Do [ M(§ + k, § - k) 
Jkn ^ ^ 

xM(-|-k,q)xGo(| + k,| + l]) 



xGo(-|-k,-|-(])Go(|-k,|-fi) 



(Al) 



where Go(q! i-^) = {—i'^ + ^i^Z^ + ^2^^^ ^ is the unrenormalized propagator. It is easy to 
evaluate the f2 integral first in flAip . Separating the f2- integral 



Go(| + k,| + fi) 



Go(|-k,|-f])rfO 



After substituting the expressions for the unrenormalized propagator in flA2l) 

z(f-l])+a 



[(f +fi)2 + 62] X [(|-fi)2 + a5 



where 



(A2) 



(A3) 



a=[^(|-M' 



6=[Ai(| + A;.)2 + Ai(| + A;, 
After integrating /n(k) over for lo — > 0, we see that, 

/n(k) 



6(a + 6) ■ 



(A4) 



(A5) 
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Substituting this fl integral in the calculation of the self-energy (1A1[) 



(2vr) 

M(-|-k,q)x ^ 



M(f + k,f-k) 



(A6) 



where a and b are defined in flA4p . Since we are interested in long- wavelength properties, 
we can do small and qy expansions. For calculating S(q, cu), we need to perform the k 



integral. Defining small parameters x = |^ and z 
in X and z 

„2 



and expanding up to lowest order 



1 



b{a + b) 2kta^ 



1 - —Ai - — ^2 tan^ 6 
2a 2a 



X 



X 



-Ai A2tan^6' + —A^ 



a 



a 



+ ^Al tan^ e+'^AiA2 tan^ 9 
a'^ a^ 



(AT) 



where Q = tan ^{^) and a = (Ai + y42tan^6'). The next step for the calculation of the 
integral is the product of two propagators M x M in (]A6I) . 

M(| + k,|-k)xM(-|-k,q) 



xzGi{Xi, A2) + (x + z)G2{\i, A2) 



+ 2G'3(Ai,A2) 

From (lATl) and (lA8l) integrand of flA6l) is, 

M(f + k,| -k) X M(-| -k,q) 



(AJ 



4 X 2A;aa2 



6(a + 6) 



X 



a 



— ^2 tan^ 9 - —Ai - —A2 tan^ 9 



a 



a 



a 



^2 

— Ai A2 tan^ 6' 

2a 2a 



X 



a 



-A, 



A2 tan^ 6' 



l2 I /i2 j.„„4 " 



+ — + — tan* ^ + ^ A1A2 tan' ^ 



a^ 



a^ 



(A9) 
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On integration (inside the [ ]) only term of O(x^), of 0(2;^) and 0(1) survive. Hence 
terms which will contribute to the integration are 



X 



a 



A, 



a 



-A2 tan^ 6 



X 



+ 2G3 1 - — - —A2 tan" 9 + —A{ 



2a 



+ —Al tan"^ 9 



2a 



a 



(AlO) 



where G2 = (2A2 + X\ — 3A1A2) and G^ = (A2 — A1A2). = kcos9 and = ksin9 and 
a = {Al + 742tan^6'). After performing the integration for these two types of terms in 
(jAini), 



E(q,aj^O) 



47r 



G'2(Ai, A2) 



{Aiql + A2ql 



+ 



G3(Ai,A2)AiA2 



This is the expression for the self-energy as given in (l29ll . 



(All) 



APPENDIX B: VERTEX RENORMALIZATION 



Here we calculate the three-point symmetrised vertex function F. There are three distinct 
one- loop diagrams r^, T^y and T^. contributing to the correction to the vertex as shown in 
(Fig[2](b)). These diagrams all have multiplicity 4. In this Appendix we will go into the 
details of the calculation of Fq. The calculations for F;, and Fc are the same as for F^. 
Small variables x and defined in Appendix \^ for self-energy. We start from the 

symmetrised version of 



F,(q,ki)=4x2Z}o / M(§ + k, § - k) 
xM(| + ki,k-ki) 
xM(|-ki,-k + ki) 



X 



Go(k-ki,fi-fii; 



Go(| + k,| + ^])xGo(|-k,|-^]) 



(Bl) 
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Separating the Q integral part from the full integration in (IBip 



+00 



Go{k-ki,Q-ni] 



X 



Go(| + k, I + n)Go{^ - I - ^)^^ 
for uj — > and Qi — > limit and writing in terms of real and imaginary parts, 



(B2) 



(B3) 



For UJ 
and 



(fi2 + 52)(^2 + a2)(fi2 + c2) 

and Qi — > limits /m(/^(k)) = 0. where a and b defined in (lA4l) 



Performing the integral over ^2, 



7i{2c + a + b) 
c{a + c){b + c)(a + b) 



(B4) 



(B5) 



Similarly for Tf, and r^, 



TT 



/rn(k) 



a(a + c)(a + 6) 

TT 



6(6 + c)(a + 6) 

Substituting this /^(k) from (IBSP in the calculation of Ta, 

r,(q,ki)=4x2Do / M(^ + k,§-k) 



(B6) 



q 



xM(| + ki,k-kO 
xM(|-ki,-k + ki: 



X 



7r(2c + a + 6) 



(Yi7) 

c{a + c){b + c){a + b) 

We are interested in long wavelength properties. By defining the small quantities x = 

k 

— and Zi = where = kcosO and k^ = ksinO, up to lowest order in 



3± 

fez' 



Xi 



X, z, xi and zi, 



7i{2c + a + b) 
c{a + c){b + c){a + b) 

7T 



l + ^Ai + ^Astan^^ 
a a 



+ ^AiA2 tan^ 6* + ^^^^AiAg tan^ 6^ 



2^2 



:b81 
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The next step for the calculation of the integral is the product of three propagators 
M X M X M 

M(| + k,|-k) xM(| + ki,k-ki) 
xM(|-ki,-k + ki) 



2(— - xizi) 



+ ('^1 -14x1.1+4x1 + 4.1 

^2 \ ^ / 3xz 



+ 1 — 1 1 ~ + S^i^i - 2a;i -2.1 + 1 
From (IBSP and (]B9p . the product inside the integral for ra(q, ki) is 

"(2^ + " + ^) xM(^^+k,^-k) 



c(a + c)(6 + c)(a + 6) 2 
M(| + ki, k - ki) X M(| - ki, -k + ki^ 

-^A^tan^^ 

a. 

H —Ai tan^ 6^ ^^1^2 tan^ d 



14xi.i 



AiA2tan^6' 



a2 



A2\ / 3x2 6x1.1 . 

— — — + 6xizi Ax 

2 y V 4 a 

A2 tan^ 6* H ^^1^2 tan^ Q 



(B9) 



+ li^AiA2tan2^^ (BIO) 

We display only those terms which give a nonzero contribution after integrating over k. 
Similarly we can obtain expressions for and 
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The total F = + Ff, + Fc = Fq + 2Ff,. After doing the integration over k, the final 
expression for F, 



F(q,ki) =2(A;,.A;, 



+ 



Stt 



Ai 
2 

2 



A2 
2 

Ai 
2 



xz 
~ 16^ 
xz 
32^ 
7xz xiZi 
32^ ^ "StT 



(Bll) 



The bare vertex is 



Fo(q,ki) = 2{k^k^ 
Ao / xz 



\\ ( xz 



(B12) 



f - xi2;i) and + Xi2;i), we 



Decomposing expression in ( ]B12[) into parts of the form 
get the corrections to ^ and Hence with this decomposition (lB12p can be rewritten 

as 



F(q, ki) = 1{kjz^) ( — - XxZx 



xz 



A1A2 



4 X 2 X Stt 



SAnAl 



2x8x87r 2x8x87r 



xz 



+ 2{k^k^) y— + xizi 
4A2Ai , 



6A^ 



(B13) 



2x8x87r 2x8x87r 
Comparing with the expression for the original vertex, the corrections to ^ and ^ are 



Ai = Ai 
A2 = A2 



\ Fi(Ai,A2)/ 

2 X 87r 
, i^2(Ai,A2)r 



2 X 87r 

where functions -Fi(Ai, A2) and -F2(Ai, A2) are defined in ([36 



(B14) 
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APPENDIX C: NOISE STRENGTH RENORMALIZATION 



Here we will compute the leading-order correction to the noise strength. The relevant 
diagram which will contribute to the integral is shown in (Fig[2|^c)); it has multiplicity of 
2. Calculating the integral with this symmetrised vertex, 



AD = 2x {2D(i 



q 



q 



M(f + k,f-k) 



Go(| + k,| + fi) 



M(-- - k,k 
^ 2 ' 2 

Go(|-k,|-fi) 



Separating the ^2 integral from the full integration and taking u 

71 



0, 



(Cl) 



(C2) 



ab{a + b) 

Expanding as in the calculation of the propagator in terms of small variables x and 

z, the terms which will contribute to lowest order are of order 1. Hence to lowest order, 



(C3) 



ab{a + b) 2k^a^ 

The next step of the calculation of the integral is the product of two propagators, M x M. 
To lowest order. 



Mq + k, I - k)M(-| - k, k - 1) = kikiix^ - \,r 



(C4) 



The final expression for the product 

^ X + k, ^ - k)M(-^ - k, k - ^1 

ab{a + b) ^2 '2 ^^2 ' 2' 

2kla^ 

After performing the integration over k in the integral fICip . 



This gives 



AD 



D = D, 



87r(AiA2)3/2 
2 X Stt 



(C5) 



(C6) 



(C7) 
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APPENDIX D: AN INCOMPRESSIBLE ACTIVE NEMATIC 



In this section we give the equation for the angle field 6, obtained from an incompressible 
velocity field v (V • v = 0). From ([2]), imposing p = constt and V ■ v = 0, and defining 
the transverse projector P = (1 — qq), we see that 



-r-ip- (V- Q) 



(Dl) 



writing Q in terms of 6 



v = -T-'p-{d,e,d,e) (D2) 

Substituting the expression for v in ([3]) to linear order in 6 the equation of motion 

G'o^(q,cj)6'q,^ = fe{q.,uj) - / 6'k,n6'q-k,c^-n I 7i 

Jk,n LV 



-^[P22(k)+P22(q-k)-Pn(k) 
-Pn(q-k)]) X (^M(k,q-k)^ 

+ 72 APi2(k) + Pi2(q - k)]k ■ (q - k) 



+ [Pn(k) + P22(q-k)]A:,(g,-fc,; 
+ [P22ik) + Pn{ci-k)%{qy-ky] 



(D3) 



where fg{ci,uj) is Gaussian random nonconserving noise with noise-noise correlation as 
defined in f l25|) . (^(^"^(q, u;) is inverse propagator, defined by 



. , «0 2 
tUJ + yq 



+ AiPn(q)g,'-A2P22(q)g' 



(D4) 



M(k, q — k) as defined in I^M), -Pulq), -P22(q) are diagonal components and Pi2(q) is the 
off-diagonal component of projection operator. We have not studied further the properties 
of this equation. 
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